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ABSTRACT 

We present the first observations of a Herbig Ae star with a circumstellar disk by the Far Ultraviolet 
Spectroscopic Explorer (FUSE), as well as a simultaneous observation of the star obtained with the 
Hubble Space Telescope Space Telescope Imaging Spectrograph (STIS). The spectra of AB Aurigae show 
emission and absorption features arising from gasses that have a wide range in temperature, from hot 
O VI emission to cold H2 and CO absorption. Emissions from the highly ionized species O VI and 
C III present in the FUSE spectrum are redshifted, while absorption features arising from low-ionization 
species like I, N I, and Si II are blueshifted and show characteristic stellar wind line-profiles. We find 
the total column density of H 2 toward AB Aur from the FUSE spectrum, V(H 2 ) = (6.8 ± 0.5) x 10 19 
cm~ 2 . The gas kinetic temperature of the H 2 derived from the ratio N(J = l)/N(J = 0) is 65 ± 4 K. 
The column density of the CO observed in the STIS spectrum is V(CO) = (7.1 ± 0.5) x 10 13 cm~ 2 , 
giving a CO/H 2 ratio of (1.04±0.11) x 10~ 6 . We also use the STIS spectrum to find the column density 
of H I, permitting us to calculate the total column density of hydrogen atoms, the fractional abundance 
of H 2 , and the gas-to-dust ratio. 

Subject headings: stars: individual (AB Aurigae) — stars: atmospheres — circumstellar matter 



1. INTRODUCTION 

AB Aurigae (HD 31293) is one of the brightest and best 
studied Herbig Ae stars, which are considered to be pre- 
main-sequence stars of intermediate mass (about 2-10 
M Q ) (Waters & Waelkens 1998). The star (spectral type 
AOVe+sh) is about 2 Myr old and is located at a distance 
of 144l 2 7 pc from the Sun (van den Ancker et al. 1998). 
Like most Herbig Ae/Be stars, AB Aur is surrounded by 
circumstellar gas and dust, some of which is distributed 
in a disk about the star. Standard stellar theory predicts 
that Herbig Ae/Be stars should not have convective lay- 
ers, and therefore, should not have chromospheres, coro- 
nae, or strong stellar winds. However, AB Aur and most 
other Herbig Ae/Be stars do show chromospheric and wind 
features (Bohm & Catala 1993). The Far Ultraviolet Spec- 
troscopic Explorer (FUSE) provides high-resolution access 
to a portion of the AB Aur spectrum previously observed 
only at low spectral resolution with the Hopkins Ultraviolet 
Telescope. In this Letter, we present FUSE observations 
of AB Aur, including the first observation of hot O VI gas 
(T ks 3 x 10 5 K) and cold molecular hydrogen. The O VI 
emission should provide an important constraint on any 
models of the stellar activity seen in the AB Aur system. 
We also present a simultaneous observation of AB Aur 
with the Hubble Space Telescope (HST) Space Telescope 
Imaging Spectrograph (STIS); this observation strongly 
complements the FUSE observations by providing column 



densities of H I and CO. 

2. fuse OBSERVATIONS AND DATA REDUCTION 

AB Aurigae was observed with FUSE on 2000 February 
27 and 28, for a total exposure time of 13.7 ks. The data 
were obtained using the low-resolution aperture (LWRS), 
and alignment of the four channels (LiF 1, LiF 2, SiC 1, 
and SiC 2) was maintained throughout the whole observa- 
tion (see Sahnow et al. (2000) for a more detailed expla- 
nation). The data were flux calibrated using calf use 1.6.8 
pipeline processing software, but more recent wavelength 
solutions (2000 June 19) were used for our analysis. We 
obtained a S/N of about 5 per 12-pixel resolution element 
in the LiF la segment (near 1060 A), and about 7 in the 
LiF lb segment (near 1140 A). 

The wavelength solutions used provide good relative cal- 
ibration across the LiF channels. However, they have sig- 
nificant zero-point offsets. In order to establish the abso- 
lute wavelength calibration of our spectra, we made use of 
the simultaneous STIS observation of AB Aur described in 
§5. Assuming that our H 2 gas is located in the same region 
and is at the same velocity as the carbon monoxide gas ob- 
served in the STIS spectrum, we set vn 2 — vco = 10 ± 3 
km s^ 1 , which is close to the heliocentric radial velocity of 
the star (v ra d = 8 km s^ 1 ). 

3. ANALYSIS OF fuse SPECTRUM 
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The spectrum in Figure la shows signal from AB Aur 
down to C III at 977 A, and is rich in emission and absorp- 
tion features, many of which we have not yet been able to 
identify. The stellar photosphere flux drops off sharply at 
around 1270 A (Figure lb), so virtually all the observed 
flux in the FUSE bandpass is due to a forest of emission 
lines and possibly some far-UV excess continuum of un- 
certain origin. 

3.1. Emission from highly-ionized species O VI and C III 

The O VI emission doublet in the LiF la segment spec- 
trum is shown in Figure 2a. The A1038 line is suppressed 
by H2 absorption; two narrow H2 lines are also superim- 
posed on the A1032 line. These two narrow lines were 
masked out of the data, and a Gaussian model was least- 
squares fitted to the spectrum between 1029 A and 1035 
A. A model for the O VI doublet was created using this 
Gaussian fit, assuming the optically thin line ratio 2:1 to 
model the A1038 line. This O VI model was multiplied 
by the normalized H 2 absorption model discussed in §3.2; 
the result is overplotted on the data in Figure 2a, show- 
ing how the strong H2 absorption lines near 1037 A and 
1039 A cut off the wings of the A1038 line, reproducing the 
peculiar line shape quite well. However, this analysis indi- 
cates that there should be additional flux visible at 1038.8 
A; absorption from some other atomic species may be re- 
sponsible for the lack of observed emission. The Gaussian 
fit to the A1032 line is redshifted (v = 105 ± 11 km s" 1 , 
corrected for the zero-point offset of the wavelength cali- 
bration), with a FWHM = 356 ± 11 km s" 1 . 

The O VI line profile does not appear to be a standard 
type I P Cygni profile (Beals 1950), because no blueshiftcd 
absorption is observed on the small but significant contin- 
uum flux around the emission line. While the Gaussian 
model fits the red wing of the A1032 line quite well, there 
is excess emission on the blue wing. This asymmetry in 
the line-profile suggests that the O VI emission is formed 
in AB Aur's stellar wind, but this conclusion is tentative 
until more detailed modeling is complete. 

The C III A977 line in the SiC lb segment spectrum is 
shown in Figure 2b. Again, there are narrow H 2 absorption 
lines superimposed on the emission feature, which were 
masked out before least-squares fitting a Gaussian to the 
line. The bulk of the emission is redshifted (v — 192 ± 11 
km s _1 ), but the C III emission feature clearly is not well- 
described by a single Gaussian. There is an "absorption 
feature" at ~ 978 A, superimposed on the C III emission, 
that does not arise from H 2 . The origin of this feature 
has not yet been determined, but it is not likely to be 
interstellar, because of its large redshift and width. 

3.2. Molecular hydrogen 

The H 2 absorption features visible in our FUSE spec- 
trum do not show any evidence of multiple velocity compo- 
nents. The LiF la and LiF 2a data were normalized with 
simple continua found from linear fits to the spectrum be- 
tween 1040 A and 1120 A. Our H 2 absorption model was 
generated using molecular data from Abgrall ct al. (1993a) 
and Abgrall et al. (1993b). Voigt line-profiles were used to 
create transmission functions, which were then convolved 
with a Gaussian instrumental line-spread function with 
FWHM = A/15000. The X 2 statistic between the model 



and the data was minimized to find the best model param- 
eters, with 1 — a error bars determined from the contours of 
X 2 . The velocity of the gas was found by least-squares fit- 
ting Gaussians to a number of narrow lines between 1050 
A and 1120 A. We measured w H2 = 103 ± 10 km s _1 , 
but set this value to the STIS carbon monoxide velocity, 
v co = 10 ± 3 km s _1 , to determine the zero-point offset 
of our wavelength calibration. 

Since the relative populations of H 2 energy levels with 
v = and J = 0, 1 are determined primarily by thermal 
collisions, the excitation temperature found from these lev- 
els should be close to the kinetic temperature of the gas. 
However, the populations of the higher levels (J > 2) can 
be inflated by UV and/or formation pumping and by ra- 
diative cascade (Shull & Beckwith 1982). Therefore, the 
low J and high J levels were analyzed separately. Win- 
dows containing only lines arising from the J = 0,1,2 
levels in the five bands between 1048 A and 1120 A were 
cut out of the normalized LiF la and LiF 2a data and 
X 2 minimization performed using a model containing only 
these lines. We permitted the model to have a slightly 
different velocity shift for each of the windows, to elim- 
inate any systematic errors arising from the uncertainty 
in the H 2 velocity determination. The parameters of the 
model were N(J = 0), N(J = 1), N(J = 2), the col- 
umn densities in the J = 0,1,2 levels respectively, and 
b, the Dopplcr broadening parameter. From the ratio 
N(J = 1)/N(J = 0), we calculate the gas kinetic tem- 
perature of the H 2 , T = 65 ± 4 K. To analyze the high J 
lines, we again cut windows containing only lines arising 
from J > 2 out of the data, and performed x 2 minimiza- 
tion using a model containing only those lines. We find the 
excitation temperature of the high J lines to be 212 ± 20 
K. The column densities of molecules in each rotational 
level were summed to find the total H 2 column density, 
AT(H 2 ) = (6.8 ± 0.5) x 10 19 cm" 2 . The final best two- 
temperature model is shown overplotted on a portion of 
the LiF la spectrum in Figure 3a. 

4. STIS OBSERVATION AND DATA REDUCTION 

AB Aur was observed by HST on 2000 February 28. 
A total exposure time of 1813 s was obtained using the 
E140M cchclle grating and the 0."2 x 0."06 spectroscopic 
slit. The spectrum covered 1150-1725 A with a resolution 
of w 46,000. The data were reduced using the STIS IDT 
software calstis. Echcllc spectra are subject to interordcr 
scattered light; to compensate for this, we used the algo- 
rithm developed by Don Lindler and discussed in Holberg 
et al. (1999), which accurately achieves zero flux in the 
trough of the interstellar Lya absorption. The STIS wave- 
length scale in this mode is accurate to 3 km s _1 (Pagano 
et al. 2000). The spectrum between 1200 A and 1300 A 
is shown in Figure lb. Longward of 1300 A, the spectrum 
resembles an early A star with the heavy line blanketing 
typical of Herbig Ae stars in the UV. The STIS spectrum 
reached a continuum S/N=20.6 near 1471 A, and S/N=16 
near 1504 A. The S/N at shorter wavelengths is lower, S/N 
< 3.8 per 0.2 A bin at Lyman a. 

5. ANALYSIS OF STIS SPECTRUM 

The STIS spectrum of AB Aur is rich in emission lines 
from a wide range of ionization states. The data short- 
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ward of 1310 A show type I P Cygni profiles arising from 
the low-ionization species O I and N I, and more complex 
emission profiles in Si III and C III A1176. 

5.1. Extinction, Lyman a, and N (H I) 

The measured color index of AB Aur, (B-V), is 0.12. 
IUE low-dispersion spectra of AB Aur (SWP 24389) 
were unreddened using the interstellar excinction law of 
Cardelli, Clayton, & Mathis (1989). Comparison with 
unreddened spectral standards suggests that the best over- 
all match in the FUV is 7 UMa, type AO V (SWP 8198), 
with (B-V)o = 0.04 and R = 3.1. We therefore find the 
selective extinction toward AB Aur, E(B — V) = 0.08, and 
the general extinction, Ay = R x E(B — V) = 0.25. 

Lyman a emission is formed in the chromosphere of AB 
Aur, and it appears likely that the chromosphere is ex- 
panding, based on profiles of other low-ionization species. 
Since these line-profiles in the STIS spectrum of AB Aur 
are observed to be type I P Cygni profiles, with the emis- 
sion extending down to the stellar radial velocity, it is 
likely that the intrinsic Lyman a profile is a type I P Cygni 
profile as well. However, the short wavelength edge of the 
observed emission does not extend down to the stellar ra- 
dial velocity; this is presumably due to the interstellar 
Lyman a absorption. The wavelength at which the inter- 
stellar absorption is no longer optically thick, where the 
observed flux rises above the zero level, is highly sensitive 
to the interstellar H I column density, but insensitive to 
the assumed b-value. We find N(R I) toward the star to 
be 3 x 10 20 cm" 2 . We obtain roughly the same value using 
several different assumptions about the underlying contin- 
uum flux. However, uncertainty about the exact shape of 
the stellar line-profile could cause the error in this value 
to be as much as 50%. Future work will refine N(R I) 
after detailed modeling of AB Aur's winds determines the 
unattenuated stellar flux presented to the interstellar H I 
atoms. 

5.2. Carbon monoxide 

The (1-0), (2-0), and (3-0) bands of the CO Fourth Pos- 
itive band system (A *n - X X E+) are visible in the STIS 
spectrum (Figure 3b). No lines arising from 13 CO were 
observed. The continua in the vicinity of the bands were 
fitted with fifth-degree polynomials. An absorption model 
was generated using wavelengths and oscillator strengths 
from Morton & Noreau (1994) and energies of the ground- 
state levels were calculated using the Dunham coefficients 
from Farrenqet al. (1991). x 2 minimization was performed 
on all three bands simultaneously to find the best rota- 
tional excitation temperature, Tco — 7.0 ±0.7 K, Dopplcr 
broadening parameter, b — 3.8 ± 0.4 km s _1 , column den- 
sity of 12 CO, AT(CO) = (7.1 ± 0.5) x 10 13 cm" 2 , and ve- 
locity centroid, v = 10 ± 3 km s _1 . 

6. DISCUSSION 

We find a total column density of hydrogen atoms, 
N H = N(R I) + 2N(R 2 ) = 4.4 x 10 20 cm" 2 , which is typi- 
cal of a diffuse molecular cloud. The fractional abundance 



of H 2 , / = 27V(H 2 )/N H = 0.31. Only 11 out of 91 stars 
in the Savage et al. (1977) Copernicus survey of interstel- 
lar H2 show a larger value of /, and these sight-lines all 
have larger values of E(B — V) and Nh. This might sug- 
gest that our linc-of-sight to AB Aur has an anomalously 
high molecular fraction, but this conclusion is tentative 
until further FUSE H 2 measurements of a variety of in- 
terstellar environments are complete. The CO/H2 ratio 
is (1.04 ± 0.11) x 10~ 6 . We find the gas-to-dust ratio, 
Ah/E(B - V) = 5.5 x 10 21 . The mean ratio for standard 
clouds within 2 kpc of the Sun is 5.9 x 10 21 , indicating that 
our gas-to-dust ratio is not atypical (Bohlin et al. 1978). 

Optical coronagraphic images of AB Aur show an ex- 
tended arc of reflection nebulosity to the east of the sys- 
tem, in addition to the circumstellar disk centered on the 
star (Nakajima & Golimowski 1995; Grady et al. 1999). 
Nakajima & Golimowski (1995) suggest that such nebu- 
losities are the remnants of the original molecular clouds 
that collapsed to form low-mass stars. Our values of the 
various parameters describing the cold H2 and CO gasses 
(Tk in , 7V H , /, CO/H2 and 7V H /E(B - V)) are consistent 
with this suggestion, and we believe that the H 2 and CO 
gasses toward AB Aur are probably located in a remnant 
molecular cloud envelope around the star. The fact that 
the velocity of the CO gas is equal (within measurement 
error) to the stellar radial velocity supports this conclu- 
sion. The molecular gasses are not likely to be located in 
the circumstellar disk around AB Aur, since the disk in- 
clination is such that our line-of-sight to the star does not 
pass through the disk (Grady et al. 1999). 

Turning to the hot, highly-ionized gasses, FUV emission 
lines from O VI and C III are observed in the solar spec- 
trum, and are produced in the transition zone between the 
chromosphere and the corona. However, detailed model- 
ing of the N V (T « 1.4 x 10 5 K) emission previously 
observed from AB Aur with GHRS found that a homoge- 
neous, spherically-symmetric high-temperature zone could 
not produce the observed N V emission without producing 
strong N IV or C IV emission as well, which is not seen 
(Bouret et al. 1997). Bouret et al. (1997) concluded that 
small, hot clumps, formed by shocks in a stellar wind, 
could reproduce the FUV emission lines, as well as the 
weak X-ray flux observed from AB Aur. Future work will 
determine if the wind-shock model can reproduce the ob- 
served emission from C III and the higher-temperature 
(T « 3 x 10 5 K) species O VI. 

This work is based on observations made with the 
NASA-CNES-CSA Far Ultraviolet Spectroscopic Explorer 
and on NASA-ESA Hubble Space Telescope observations 
obtained at the STScI under the Guaranteed Time Ob- 
server program 8065. FUSE is operated for NASA by the 
Johns Hopkins University under NASA contract NAS5- 
32985. The work performed by J.-C. Bouret was supported 
by a National Research Council-(NASA GSFC) Research 
Associateship. 
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Fig. 1. — Overview of the AB Aurigae FUV spectrum. The 1 — cr flux measurement errors are plotted with a dashed line. Airglow lines 
are indicated by ©, and prominent features arc labeled, (a) FUSE spectrum of AB Aur, combining data from all channels. The data have 
been rcbinncd by a factor of 80 for this plot. (b) Short wavelength portion of the STIS E140M spectrum of AB Aur. The data have been 
rcbinned by a factor of 10 for this plot. 
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Fig. 2. — Emission from highly-ionized species in the FUSE spectrum. The data have been rebinned by a factor of 8 for these plots. Our 
best H2 absorption model is shown at the top of both plots. The 1 — a flux measurement errors are plotted with a dashed line. (a) O VI 
AA1032, 1038 doublet in the LiF la spectrum. An optically thin Gaussian O VI model is shown with a dotted line. The normalized H2 model 
was multiplied by the O VI model and overplotted with a thick line. (b) C III A977 line in the SiC lb spectrum. Our Gaussian fit to the 
line multiplied by the H2 model is overplotted with a thick line. 
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Fig. 3. — Molecular gasses toward AB Aur. (a) H2 Lyman (3-0) band in FUSE LiF la spectrum. The data have been rebinned by a 
factor of 10 for this plot. The 1 — a flux measurement errors are plotted with a dashed line. The best two-temperature model is ovcrplotted 
with a thick line, (b) Three bands of the CO Fourth Positive band system in the STIS E140M spectrum. The data have not been rebinned 
or smoothed for this plot. The 2 — <r flux error bars are overlaid on the data. The best model is overplotted with a thick line. 



